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Strongly correlated amorphous solids are a class of glass-formers whose inter-particle potential 
admits an approximate inverse power-law form in a relevant range of inter-particle distances. We 
study the steady-state plastic flow of such systems, firstly in the athermal, quasi-static limit, and 
secondly at finite temperatures and strain rates. In all cases we demonstrate the usefulness of 
scaling concepts to reduce the data to universal scaling functions where the scaling exponents are 
determined a-priori from the inter-particle potential. In particular we show that the steady plastic 
fiow at finite temperatures with efficient heat extraction is uniquely characterized by two scaled 
variables; equivalently, the steady state displays an equation of state that relates one scaled variable 
to the other two. We discuss the range of applicability of the scaling theory, and the connection to 
density scaling in supercooled liquid dynamics. We explain that the description of transient states 
calls for additional state variables whose identity is still far from obvious. 



I. INTRODUCTION 

The equations of fluid mechanics appear to provide an 
adequate description for the flow of liquids for an ex- 
tremely vifide range of boundary conditions and exter- 
nal forcing. A similarly successful theory is still lacking 
for the description of elasto-plastic dynamics in amor- 
phous solids which form as the result of the glass tran- 
sition. While being essentially "frozen liquids" , amor- 
phous solids differ from regular liquids in having a yield 
strength as , a material parameter which depends on the 
density, temperature etc, which is the maximal value of 
the internal stress that the material can support by elas- 
tic forces. Regular liquids cannot support any amount 

\ of stress without flowing. When the stress exceeds the 
yield strength the material begins to respond plastically, 
and under a given external shear rate can develop a 
steady state plastic flow with a mean "flow stress" <Joo- 
The analog of the Navier-Stokes equations which can de- 
scribe the whole spectrum of elasto-plastic responses in 
terms of macroscopic variables is not known yet, and 

\ their derivation is the subject of much current research 
[l, 2, 3, 4, 5, with significant amount of debate. 

In this paper we focus attention on the steady-state plas- 
tic flow which is obtained under the action of a constant 
external strain rate. We will argue below that the char- 
acterization of such a state is considerably simpler than 

' the full description of transient states, the latter call for 
a larger number of macroscopic variables whose nature is 
not obvious and the constitutive relations between them 
are not known. For the steady plastic flow state we can 
make progress and determine what are the state variables 
that determine the state uniquely. 

To simplify things further we limit our attention at 
present to materials whose inter-particle potential can be 
approximated, for the range of inter-particle distances 
of relevance, by an inverse power law potential. This 
same class of materials and the interesting scaling prop- 
erties that they exhibit attracted considerable interest 



in the context of the dynamics of super-cooled liquids, 
first experimentally [M Q [S [3 

and then theoreti- 
cally [13, m, [H, HO] ■ In the context of the mechanical 
properties of amorphous solids we believe that the first 
example of using the special scaling properties of these 
materials appeared in Q where focus was put on the 
athermal limit and quasi-static strain. In this paper we 
explore further the quasi-static limit, and then extend 
the discussion to systems at finite temperatures and fi- 
nite strain rates. The discussion culminates with finding 
which are the minimal number of re-scaled state variables 
that determine uniquely the steady plastic flow in such 
materials. Any general theory that attempts to provide 
a complete description of elasto-plasticity in amorphous 
solids should reduce, in the steady flow state of materials 
of the present class, to a theory that contains these and 
only these variables. 



The structure of the paper is as follows: In Sect. [Ill we 
introduce the systems under study, and explain how they 
are simulated both in the athermal, quasi-static limit and 
at flnite temperatures and strain rates. In Sect. IIIII we 
explain the special scaling properties that these systems 
possess, and predict theoretically what is expected in the 
steady plastic flow state. This is the central part of the 
paper. We then provide detailed presentations of simu- 
lation results and demonstrate how they compare to the 
predictions of the scaling theory. We discuss analytic 
properties of the scaling function, and demonstrate the 
conditions under which the scaling breaks down. In Sect. 
IIVI we discuss the consequences of our thinking to super- 
cooled liquids, and propose that the scaling function used 
in the literature in this context are incomplete. Sect. |V| 
summarizes the findings, and provides a discussion of the 
road ahead, especially in terms of extensions to transient 
states. 
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II. SYSTEMS AND METHODS OF 
SIMULATION 

A. System Definitions 

In this work we employ two-dimensional polydisperse 
systems of point particles of equal mass to, interacting 
via two qualitatively different pair- wise potentials. Each 



particle i is assigned an interaction parameter Ai from a 
normal distribution with mean (A^) = 1. The variance 
is governed by the poly-dispersity parameter A — 15% 
where = ii^^^^W)!!. With the definition Ay = ^(Ai + 

Aj) the first potential UR{rij) is purely repulsive, of which 
the shape is characterized by the interger k: 
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We chose Bq = 0.2 for all systems discussed, and vary 
the integer k in the following. This pair-wise potential is 
constructed such as to minimize computation time, and 
is smooth up to second derivative, which is required for 
minimization procedures. 

The second pair- wise potential {/^(r^ ) reads 



with U {vij ) = 




r.{X,,) < r < r,(A,j) (2) 
r > rc(Ay ) 



k = 12, 



and Tc = l.SGA^j. The attractive part U{r) 
is glued smoothly to the repulsive part. We choose 

= hP {77^) ^here P{x) = Eto and the 
coefficients Ai (see Table |l| are chosen such that the po- 
tential is smooth up to second derivative. These pairwise 



attractive UA{rij), k - 
repulsive UR{rij), k = 8 
repulsive UR{rij), k = 10 




FIG. 1: Color online: The different pairwise potentials discussed 
in this work. 
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TABLE I: The coefficients in P{x) = ELo^'^'- ^ee text. 



instantaneous quenching of random, high temperature 
configurations; this explains the apparent noise and ab- 
sence of stress peaks in the transients. Furthermore, it is 
important to note that due to finite system sizes, the ini- 
tial value of the stress of the quenched configurations in 
some experiments is non-zero; this is however irrelevant 
for steady state statistics. 



potentials are displayed in Fig.[T]for the cases of interest. 
Below the units of length, energy, mass and temperature 
are A = (Ai), e, to and e/ks where ks is Boltzmann's con- 
stant. The time units r^, are accordingly t,^ = yJmX^ /e. 
From here and in the following we denote the density as 
p = Y^ ^'^d define the dimensionless density p = X^p. 
Also, we will refer to the dimensionless density as just 
the density, for the sake of brevity. 

Initial conditions for all the simulations, for both meth- 
ods described in the next Subsection, were obtained by 



B. methods 

The work presented here is based on two types of sim- 
ulational methods. The first type corresponds to the 
athermal quasi-static (AQS) limit T — > and 7 — > 0, 
where 7 is the strain rate. AOS methods have been ex- 
tensively used recently 0, 0, [a, 0j H, Q as a tool for in- 
vestigating plasticity in amorphous systems. The order 
in which the limits T — s- 0,7 — > are taken is impor- 
tant, since one expects that at any finite temperature 
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the stress in the system can thermany relax given long 
enough time (or small enough strain rates), hence 
the limit T ^ should be taken prior to the 7 — > 
limit. According to AQS methods, starting from a com- 
pletely quenched configuration of the system, we apply 
an afhnc simple shear transformation to each particle i 
in our shear cell, according to 

fiy ^ Tiy j (3) 

in addition to imposing Lees-Edwards boundary condi- 
tions The strain increment 5e plays a role analo- 
gous to the integration step in standard MD simulations. 
We choose for the discussed systems 5e — 10^"*, which 
while not sufhciently small for extracting exact statistics 
of plastic flow events as done in Q, it is, however, suf- 
flciently small for the analysis of the steady state prop- 
erties and mean values. The affine transformation ([3|) 
is then followed by the minimization ^12i] of the poten- 
tial energy under the constraints imposed by the strain 
increment and the periodic boundary conditions. We 
chose the termination threshold of the minimizations to 
be |V?7|V7V = 10-1^ 

The second simulation method employs the so-called 
SLLOD equations of motion [llj]. For our constant strain 
rate 2D systems, they read 

Tix = Pix/m + , 

^iy — Piy I ^ i 
Pix — fix '~iPiy ; 
Piy = fiy ■ 

We use a leapfrog integration scheme for the above equa- 
tions, and keep the temperature constant by employing 
the Berendsen thermostat [ll|, measuring the instanta- 
neous temperature with respect to a homogeneous shear 
flow. The integration time steps were varied between 
5t = 0.007 and 5t = 0.001, depending on density, such 
that numerical stability was maintained for all densities 
simulated. The time scale tt for heat extraction [llj was 
chosen such that rate of heat generation is smaller than 
the rate of heat extraction. For the lowest densities this 
was chosen to be tt ~ lOr^,. 

III. THE SCALING THEORY 




FIG. 2: Color online: r-^^^ in the range of ro/A G 

[pmix^, p~i/^] for A: = 8 in green asterisks, and for = 10 
in blue circles. The line through the points represents the 
scaling laws (O. 

the range of densities employed here is equivalent to a 

range of tq/A G [pmix'^, Pmin ''l • We flnd that in this range, 
to a very good approximation, 

r'^-i dr ^ Xd \\) ■ ^ > 

In two dimensions v — 4.80 for fc = 8 and v ~ 5.87 for 
k = 10, see Fig. H 

In the following discussion we deflne the flow stress 
(Too to be the steady-state value of the stress under con- 
stant external strain rate. In general, the flow stress is 
a function of a set of state variables, which specify the 
conditions in which the experiments are carried out. For 
the systems and experiments discussed in this work, the 
flow stress depends on the density p, the temperature 
T, and the strain rate 7. In addition, one can expect 
also a dependence on the heat extraction rate . We 
choose to exclude the latter from the present discussion, 
and we do so by choosing the rate of heat extraction to 
be much larger than the rate of heat production. So, we 
propose at this point that a 00 — <Joo {T, p, 7) . The yield 
stress (Ty{p) is defined as the steady state value of the 
stress under the limits T — > and 7^0 (see discussion 
regarding these limits in Subsect. IIII Cp . i.e. 

cry = CToo(p, r ^ 0,7 0) . (6) 



The discussion of the relaxation properties of g lass for- 
mers in the super-cooled regime [l^, Q, El, 13 [13, [H, 
[Tgl . l20j and of the mechanical properties of the amor- 
phous solids 's'l simplifies significantly when the inter- 
particle potential assumes an effective inverse power-law 
from in the relevant range of inter-particle distances. As 
an example consider the potential ([1]) in the density range 
p € [1, 1.6]. Since in d dimensions the characteristic inter- 
particle distance scales like 



A. Scaling in the Athermal, Quasi-static limit 

In the athermal, quasi-static limit the only parame- 
ter left is the density; consideration of the temperature 
and strain rate effects will be taken up in the next Sub- 
section. Denote the distribution of inter-particle dis- 
tances as p{r); then the mean inter-particle distance is 
ro(/9) = / rp(r; p)dr. Note that this probability distri- 
bution only accounts for distances which are relevant in 



terms of the interaction, namely for < Ay 
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FIG. 3: Color online: stress-strain curves averaged over 20 
independent runs for an athermal system with A'^ = 4096, 
= 8 (left panel) and k — 10 (right panel) as a function of 
the density, with the density increasing from bottom to top. 



If p(r) is sufficiently sharply peaked around ro, we can 
write 



dr 



OUr 



dr 



(?) 



(7) 



From here we predict that for our systems with short- 
range forces the scaling of the yield stress should be 



rn ^1 

cry ^ J\ — ~ j-^p 



V 



(8) 



In the athermal, quasi-static limit the shear modulus 
must obey the same scaling 



(9) 



These scaling laws lead to the expectation that re- 
plotting stress-strain curves in terms of re-scaled variable 
a/p'^ should result in complete data collapse. Indeed, 
our simulations vindicate this expectation. In Fig. [3] 
we present the raw stress-strain curves in the athermal, 
quasi-static limit using seven different values of the den- 
sity. For each density we simulated 20 independent runs 
of iV = 4096 particles, using the pairwise potential ([T]) 
and two choices of the integers fc = 8 and k = 10. Fig. 0] 
demonstrates the superb data collapse for the scaled vari- 
able. The insets are a direct test of the scaling laws ([8]) 
and (13. 



Scaling Theory with Temperature and External 
Strain Rate 



Once we perform measurements at finite temperatures 
and external strain rates the scaling considerations must 
incorporate temporal and energy scales. The typical free 
energy density in the steady-state plastic flow should 
scale like cry x Se where Se is the typical strain inter- 
val between plastic events, Se ay /a*. Accordingly, the 




FIG. 4: Color online: The same stress-strain curves as in 
Fig. [3] but with the stress rescaled by p" , with v = 4.80 for 
k = 8 (top panel) and v = 5.87 for fc = 10 (bottom panel). 
The insets demonstrate the density dependence of ay and p 
according to p" . 



intensive energetic contribution to barriers SG (that gov- 
ern thermal activation) scales with the density according 
to 



(SG) 





N 



(10) 



Note that this is the "density scaling" proposed in [Ij, 
Hi, [S [H, [13, [ll m, m in the context of the dynamics 
of super-cooled liquids. For the present purposes we need 
to explore further scaling relations; we estimate now the 
density scaling of the typical time-scale tq with respect 
to which all the rates in the theory should be compared. 
We begin with the speed of sound c^; using Eq. © we 
write 



We can now define the time scale tq = ra/cg 
Eqs. dU and ((Til) we obtain 

iyd-d + 2 

Tq Ti,p 2d 



(11) 

Using 
(12) 



Using Eq. [TO] we conclude that the effect of tempera- 
ture on the dynamics in the steady state must be in- 
variant once the temperature is rescaled by p'^~^. On 
the other hand the external strain rate 7 should leave 

• • vd-d+2 

the system invariant once rescaled by p 2d due to 
Eq. [T2I Putting together all these we finally propose the 
expected scaling- function form for the flow stress CTqo: 



(13) 
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This is the central theoretical result of this section. We 
stress that we chose to favor the flow stress and wrote it 
in terms of the scaling function of the other two dimen- 
sionless variables. We could equivalently choose any of 
the other two variables to be represented in an analog 
way in terms of two dimcnsionless variables. This scal- 
ing function form is in fact an equation of state for the 
steady plastic flow. 

For d ~ 2 this general result assumes the form 



T 



7 



(14) 



To demonstrate the high degree of precision with which 
the scaling theory is obeyed we performed simulations at 
flnite temperature and strain rate (see methods section) 
in which we prepared 10 independent systems (for each 
density) of iV = 10000 particles at the densities p = 
1.0,1.1,1.2,1.3 and 1.4. Defining the two dimcnsionless 



variable 



_ T 



and y = 



, we fix the value 



yo = 1.6 X 10 for all densities, and simulated all the five 
densities for the values x = 0.001,0.01,0.1 and 0.2. The 
results are displayed in Fig. ([5]). We see the excellent 




FIG. 5; Color online: Left Panels: stress normalized by p"^ 
vs. strain for the x values x = 0.001, a; = 0.01, x = 0.1 and 
X — 0.2, increasing from top to bottom. Right panels: log-log 
plots of the steady state flow stress as a function of density, 
for the same corresponding values of x. 



ulated the entire range of x values for which S{x, y) ex- 
ists. The result is shown in Fig. [HI in addition to the data 
obtained for all the other densities and x values shown 
in Fig. O The excellent data collapse is quite apparent. 
It is noteworthy that at low temperatures the function 



-AQS limit i; ^ 0, y^Q 
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FIG. 6: Color online: The function <S (a;; yo)- Data is displayed 
for p = 1.15 (blue circles) over a wide range of a; = ^^l^-i 
values, and for the densities of Fig. [5] over the x values x = 
0.001,2; = 0.01, s = 0.1 and x = 0.2. The value of yo = 
„/2 is 1-658 X 10~^ for all simulated systems. 



reaches smoothly, albeit with a very high gradient, pre- 
cisely the athermal, quasi-static limit that was studied 
in the previous Subsection. The high gradient as T — > 
in a similar, experimentally obtained function, was inter- 
preted in [2l| as resulting from quantum-mechanical ef- 
fects. Obviously in our purely classical simulations there 
are no quantum effects and it remains very interesting to 
unfathom the origin of the very fast change in the flow 
stress over a very short temperature interval. 

To emphasize the relevance of the temporal scaling we 
simulated steady fiow states at different external strain 
rates but at the same x values. The result are shown 
in Fig. [71 We see that as the temperature increases, 
the relative sensitivity of the flow stress to changes in 
the the strain rate increases appreciably. Note that the 
value of j/o = 1-6 X 10~^ for which the data collapse was 
demonstrated is well within the range of high sensitivity 
to changes in the strain rate. In other words, without 
rescaling the strain rate properly there is no hope for 
data collapse. Further analytic properties of the scaling 
function are discussed in the next Subsection. 



data collapse and also the quality of the scaling laws for 
the flow stress; the slopes of the lines in the right panels 
are those predicted theoretically in Eq. (HH), i.e. CToo ^ 

We now test the quality of the prediction of the exis- 
tence of the scaling function iS(x, y). To this aim we fixed 
a value of p = 1.15 and the same j/o = 1.6x 10~^, and sim- 



C. Analytic Properties of the Scaling Function 

The entire physics of the steady flow state for this class 
of systems is encoded in the scaling function S{x,y). It 
is therefore very challenging to derive the form of this 
functions from first principles. We are not yet in a po- 
sition to do so; at this point we can only present the 
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FIG. 7: Color online: The scaling function S{xo,y) normal- 
ized by the values S{xo, y = 2.5 x 10" 

y- 



for various values of 



analytic properties of this function as a preparation for 
future discussions. 

Firstly, it is noteworthy that the limits limj;_,o limj,_+o 
do not commute. We expect that 



and limj^^o 1™, 



lim lim y) = 

x—^O y—^Q 



(15) 



simply because at any finite temperature, given enough 
time to relax the stress, the flow stress must vanish |1C| . 
On the other hand 



lim lim S{x, y) = <7y / 



(16) 



as can be seen directly from Fig. [HI 

Secondly, in the athermal limit x — > the flow stress 
loses its dependence on the external strain rate for suffi- 
ciently small values of y, 



dS{x,y) 



lim lim 

y^Ox^a ay 







(17) 



This property can be seen directly in Fig. [T] The phys- 
ical reason for this property is that without substantial 
thermal activation the physics becomes insensitive to ex- 
ternal time scales. This limit is expected to hold when 
the external strain rate is much smaller than the elastic 
relaxation rate; interplays between high strain rates and 
the flow stress were investigated in [22]. 

Finally, we observe an inflection point in S{x,y), see 
Fig. [HI where 



dx^ 



= 



(18) 



We conjecture that this inflection point separates a "low 
temperature region" from a "high temperature region" 
in which the elasto-plastic physics is not the same. It is 
possible that a change from delocalized plastic events to 
more localized events [1, [l^l is the fundamental reason 
for this change, but further study is necessary to pinpoint 
this issue in a convincing way. 
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FIG. 8: Color online. Left panel: stress-strain curves for the 
potential ((2]) which has a repulsive and an attractive part. 
Right upper panel: demonstration of the failure of rescaling 
of the stress-strain curves. Lower panels: ay and (/i) as a 
function of the density. Note that predictability is regained 
only for higher densities, the straight line is p^. 



D. Applicability of the Scaling Theory 

At this point it is appropriate to discuss the general 
applicability of the scaling approach. It is sufficient to 
delineate this applicability in the context of the ather- 
mal, quasi-static limit using systems in which the inter- 
particle potential cannot be usefully approximated as in- 
verse power laws. In some model systems, e.g. [l9j . it 
has been shown that density scaling of the dynamics of 
super-cooled liquids still holds in spite of the presence 
of attractive forces in the potential. Furthermore, the 
same qualitative density scaling has been applied to a 
wide variety of experimental data, with substantial suc- 
cess [l^ [3 [3 • In these experimental systems there 
are definitely attractive forces between the particles, and 
thus the question of the applicability of the scaling theory 
is highly pertinent. 



1. Simulations 

We have simulated systems with the potential UA{r), 
Eq. ([2]) in the athermal, quasi-static limit. In this poten- 
tial an attractive branch is added to the repulsive one, 
see Fig. [TJ We again prepared 20 independent runs for 
each of the 7 densities p — 1.0,1.1,1.2,1.3,1.4,1.5 and 
1.6, this time for systems of iV = 2500 particles, and 
collected statistics for the steady state stress values (see 
methods section), as previously described. 

The raw data of the stress-strain curves is displayed in 
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FIG. 9: Color online: The pure number as a function of the 
density for the three potentials discussed in the text. Note 
that n appears to increase with the exponent of the repulsive 
part of the potential whenever scaling prevails. 



the left panel of Fig. [51 In the right upper panel we show 
what happens when we try to collapse the data by rescal- 
ing the stress by cry- Of course the stress-stain curves 
now all asymptote to the same value, but the curves fail 
to collapse, since (fi) does not scale in the same way as 
cry. Nevertheless, even in the present case we can have 
predictive power for high densities. When the density 
increases the repulsive part of the potential ^ becomes 
increasingly more relevant, and the inner power law r~^^ 
becomes dominant. We therefore expect that for higher 
densities scaling will be regained, and both ay and (/i) 
would depend on the density as p^. The two lower right 
panels in Fig. [5] show how well this prediction is realized 
also in the present case. 



2. Constancy of the ratio of the shear modulus and the 
yield stress 

Another way of flushing out the failure of scaling when 
there exist attractive forces is provided by the ratio 



ay 



(19) 



This is a pure number, which has been claimed to be 
universal for a family of metallic glasses 2l|. For sys- 
tems in which our scaling analysis holds, we have seen 
that the shear modulus scales with density in exactly the 
same manner as the yield stress (see Eq. (Ill),®), hence 
the number Q should be invariant to density changes, 
for a given system. However, when compared across dif- 
ferent systems, there is no a-priori reason to expect this 
number to be universal. Fig.[5]displays the measured val- 
ues of Q for our athermal, quasi-static experiments, for 
two different repulsive potentials of the form ([T]), using 
fc = 8 and k = 10, and for the attractive potential 
with k = 12. For the two repulsive potentials, we find 
from our numerics that this parameter differs by about 



5%, indicating non-universality. The lack of universality 
is even clearer with the last potential It is apparent 
that when scaling prevails the value of fi is constant up to 
numerical fluctuations. In the third case, where scaling 
fails, r2 is a strong function of p except at higher densities 
where scaling behavior is recaptured as explained. We 
can therefore conclude that the approximate constancy 
of ft found in a family of metallic glasses [2l| , is not fun- 
damental but only an indication of the similarity of the 
potentials for this family. In general f2 can depend on the 
inter-particle potential. It is quite clear from considering 
Eqs. ([7]), ([5]) and that the coefficients in the scaling 
laws (HI) and ([9]) may well depend on the exponent k in 
the repulsive part of the potential. The ratio of these 
pre-factors, being a pure number, could be independent 
of fc, and Q could be universal. It appears however that 
(/i) is increasing more with k than ay, and therefore 
shows a clear increase upon increasing k. At present this 
must remain an interesting riddle for future research. 



IV. RELATION TO DENSITY SCALING IN 
SUPERCOOLED LIQUIDS 

The destruction of scaling for low-density systems with 
the attractive potential ^ is in apparent contradiction 
to density scaling analysis of relaxation times in super- 
cooled liquids. As mentioned above, it has been shown in 
the context of the dynamics of supercooled liquids, both 
in model systems and in experiments, that the presence 
of attractive forces in the pairwise potentials can still be 
consistent with density scaling. In our context of mechan- 
ical properties scaling is regained only at high densities; 
it is desirable to understand whether there is a qualita- 
tive difference between the influence of attractive forces 
on mechanical properties, and the influence of attractive 
forces on the dynamics of supercooled liquids. 

The standard way in which density scaling is presented 
in the context of the d yna mics of supercooled liquids is 
in the form [H, HI M,M, El, [II 



(20) 



where Tq, is the a-relaxation time and T(x) is a scaling 
function of one rescaled variable; the exponent 7 corre- 
sponds to — 1 in our scaling analysis. 

In our opinion this form cannot be exact, and we pro- 
pose now an alternative form in light of the analysis pre- 
sented above. The form ((20)l account only for the density 
scaling of the free-energy barriers for thermal activation. 
We have noted above that on top of this the microscopic 
time-scale tq, with respect to which rates are compared, 
also varies with density, see Eq. and discussion in 

Subsect. ImBl 

Write the a-relaxation time in the standard transition- 
state-theory form 



,{T,p) = Tae— 



(21) 
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The free-energy barrier 6G scales with density as SG ^ 
ep"~^ (see discussion prior to Eq. ((TU]) ); the microscopic 
time scale should scale as tq ~ r*p~ ""^^.t^^ , (see discus- 
sion prior to Eq. (|12p ). Combining these considerations, 
we obtain the scaling form 

r„(r,p)-r,p . (22) 

We believe that this correct form was missed because 
the scaling of thermal activation barriers appears in the 
exponent of the RHS of (j2ip . whereas the scaling of the 
microscopic time scale is in the pre-factor. Nevertheless 
it is our suggestion that data should be re-analyzed using 
the proper form of the scaling function. 

V. SUMMARY AND THE ROAD AHEAD 

In this paper we offered some modest inroads into pro- 
viding a theory for elasto-plastic dynamics. We must ad- 
mit that a complete theory of elasto-plastic response of 
amorphous solids is still out of reach, mainly because of 
some fundamental riddles that are highly debated. Our 
proposition in this paper is that understanding the steady 
plastic flow state is firstly simpler than and secondly 
mandatory for achieving a full theory of elasto-plasticity. 
By focusing on glass formers with simple effective inverse 
power-law potentials we achieved a scaling theory for the 
steady-state flow stress under constant strain rate and 
finite temperatures. We have shown that in the ather- 
mal, quasi-static limit the yield stress exhibits power-law 
dependence on the density, as does the shear modulus. 
It was then shown that temperature and external strain 
rate can be incorporated into the scaling approach by ac- 
counting for thermal activation effects via energy scaling, 
and rate effects via temporal scaling. The finite temper- 



ature and finite strain rate theory appears in excellent 
agreement with the athermal, quasi-static limit when the 
appropriate limits are taken. 

The first task ahead is to provide an understanding 
from first principles of the scaling function S{x,y). We 
have discussed some analytical properties of this scal- 
ing function, some of which offer fascinating riddles for 
future research. Probably the most intriguing of these 
is the infiection point in S{x,y), see Eq. (IT51) and the 
corresponding discussion. Understanding the origin of 
this inflection point may shed light on the possibility of 
constructing mean field theories of plasticity at least for 
steady states, including the external parameter regimes 
for which they might be valid. 

Probably the most important remaining issue is the 
identification of additional state-variable that are neces- 
sary to describe transient states. It is well known that af- 
ter straining in one direction and reaching a steady state, 
a change in straining direction with an angle with respect 
to the original direction results in angle dependent tra- 
jectories. This means that a tensorial order parameter 
is written into the material during the steady flow state, 
and this object does not appear in our analysis. It must 
appear however in the transient trajectories. The iden- 
tification of this tensorial object will call for additional 
future work. 
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